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Abstract 

We present 10.27 hrs of photometry of the roAp star HD 134214 obtained by 
the MOST 1 satellite. The star is shown to be monoperiodic and oscillating at a 
frequency of 2948.97 ± 0.55 /iHz. This is consistent with earlier ground based 
photometric campaigns (e.g. Kreidl et al. 1994). We do not detect any of the 
additional frequencies identified in the recent spectroscopic study by Kurtz et 
al. (2006) down to an amplitude limit of 0.36 mmag (2a significance limit). 

1. Introduction 

The rapidly oscillating Ap stars (roAp) represent a unique subset of the chem- 
ically peculiar stars of the upper main sequence. In general, the Ap stars have 

1 MOST (Microvariability and Oscillations of STars) is a Canadian Space Agency 
mission, operated jointly by Dynacon, Inc., and the Universities of Toronto and British 
Columbia, with assistance from the University of Vienna. 
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globally organised magnetic fields of strengths of order kiloGauss, and spec- 
tral anomalies that are interpreted as vertical and horizontal chemical inho- 
mogeneities in the stellar atmosphere. The roAp stars (among the coolest 
members of the Ap class) exhibit rapid oscillations in photometry and spec- 
troscopy. These variations (first observed by Kurtz (1982)) have periods from 
about 5 to 20 minutes and low amplitudes (B <; 10 mmag). They are consistent 
with acoustic (p-mode) pulsations of low degree and high radial overtone. A 
thorough review of the roAp stars is provided by Kurtz & Martinez (2000). 

HD 134214 is in many ways a typical Ap star. It has an effective temperature 
of ~ 7500 K, is a moderately slow rotator with a usini « 2.0 km/s, and has a 
strong magnetic field with B r w —2.9 kG (parameters are derived in Shavrina 
et al. 2004). However, among the roAp stars, HD 134214 stands alone as 
the star with the shortest known period (~ 5.6 minutes). This high frequency 
oscillation is well above the estimated isothermal acoustic cutoff for a typical 
A-star model (see, e.g., Audard et al. 1998). There is also evidence that this 
oscillation frequency is variable at the level of a few tenths of fiHz over a time 
scale of approximately one year, even though the amplitude remained stable 
over many years (see Kreidl et al. 1994). Kurtz (1995) discuss the frequency 
variability observed in this and other roAp stars. 

Recently, Kurtz et al. (2006) have shown that in some roAp stars the spec- 
troscopic variability is dramatically different from that observed photometrically. 
Although HD 134214 is observed to be monoperiodic in photometry over many 
years (Kreidl et al. 1994), Kurtz et al. present evidence for up to 6 periodicities 
in their (admittedly short) spectroscopic time-series of this star. 

In this paper we present MOST (Microvariability & Oscillations of STars) 
spacebased photometry of HD 134214 and compare the observed oscillations 
to the results of Kurtz et al. (2006). 

2. Photometry and Frequency Analysis 

MOST is a Canadian microsatellite that was launched into a 820-km Sun- 
synchronous polar orbit in June 2003. Its primary science objective is to obtain 
nearly continuous, ultra-precise photometric measurements of stars for the pur- 
poses of asteroseismology. The instrument is a 15-cm Rumak-Maksutov tele- 
scope that illuminates a CCD photometer through a custom broadband filter 
(350-700 nm). Detailed technical information on the instrument is contained 
in Walker, Matthews et al. (2003). The first science results for the mission 
were published by Matthews et al. (2004). 

The roAp star HD 134214 was observed by MOST as a trial Direct Imaging 
target on 1 May 2006. Direct Imaging is the observing mode where a defocused 
star image is projected onto an open area of the Science CCD (see Rowe et al. 
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2006; and references therein, in this issue for more details on different observing 
modes of MOST). The star was in one of two target fields observed during 
each MOST satellite orbit (period = 101.413 min), so data are collected with 
a duty cycle of only ~ 47%. Exposures are 1.5 sec long, obtained every 10 
seconds. There were a total of 505 measurements collected during 10.27 hrs. 

Photometric reductions for Direct Imaging targets are described by Rowe et 
al. (2006b). In addition to these procedures, the data are detrended using a 
running mean of approximately 15 mins to reduce any remaining low-frequency 
contributions to stray light from scattered Earthshine which is modulated with 
the MOST satellite orbital period. This detrending does not influence the 
high-frequency domain of interest for roAp oscillations. The resulting changes 
in the mean level of the light curve are not large enough to affect the measured 
amplitude of the oscillation signal, within the uncertainty from the fit (see 
below). 

The reduced photometry is presented in Figure 1. The light curve of HD 
134214, shown in the top panel, has a standard deviation of <~ 1.8 mmag and 
a 2cr standard error of 0.2 mmag. Data subsets labelled A and B are shown in 
the lower two panels. 

Frequency determinations are made using CAPER; see Walker et al. (2005) 
and Saio et al. (2006). CAPER is an iterative procedure that identifies peri- 
odicities in Fourier space and simultaneously fits a set of sinusoids in the time 
domain. This method of determining oscillation parameters from time-series 
data follows the philosophy of the popular software packages Period98 (Sperl 
1998) and Period04 (Lenz & Breger 2005). Discrete Fourier Transforms (DFTs) 
are calculated and the fit is subtracted from the data successively, until there 
is no meaningful change in the fit residuals. 

The DFT and spectral window function for the HD 134214 data are shown 
in Figure 2. There is one periodicity identified at 2948.97 ± 0.55 /iHz with a 
signal-to-noise ratio of 15.4. The noise in the Fourier domain was estimated as 
the mean of the entire spectrum (0.12 mmag). We see no evidence for other 
periodicities at a detection threshold of about 2a (or 3x the noise; see Kuschnig 
et al. 1997). The fitted values of amplitude and phase are 1.88 ± 0.09 mmag 
and 1.97 ± 0.39 radians respectively. The uncertainties in the fit parameters 
are determined using a bootstrap technique (see below). The dark line in Figure 
1 shows the fit of this single sinusoid plotted over the observations. 

The uncertainties we calculate for time-series parameters derived from non- 
linear least-squares methods depend on the noise of the data (which may be a 
combination of instrumental and random processes) and on the time sampling. 
It is also known that fitted phase and frequency parameters are correlated, lead- 
ing to underestimated uncertainties in these parameters when calculated from 
a covariance matrix (see, e.g., Montgomery & O'Donoghue 1999). 
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Figure 1: The MOST photometry of HD 134214. The top panel shows the 10.27 hrs 
of data. Gaps in the light curve are the result of observing this target every second 
satellite orbit. The zoomed regions; labelled A and B, are shown respectively in the 
lower panels. The open circles are approximately twice the standard error in size 
(~ 0.20 mmag). The solid line represents the fit of a single sinusoid with a frequency 
of 2948.97 mHz 



The " bootstrap" is an ideal way to assess the uncertainties in fitted param- 
eters for a time-series analysis. The procedure is outlined in Wall & Jenkins 
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(2003) and has recently been used in a number of MOST applications (see, 
e.g., Rowe et al. 2006b and Saio et al. 2006). Clement et al. (1992) also used 
a bootstrap to estimate the uncertainties in Fourier parameters they derived for 
RR Lyrae stars. 

In short, the bootstrap is a technique that allows a user to produce a dis- 
tribution for each calculated parameter by constructing a large number of light 
curves from the original data. No assumptions need to be made about noise 
properties of the data and individual photometric errors are not needed. Each 
new light curve is assembled by randomly selecting N points from the original 
light curve (also containing N points) with the possibility of replacement. In this 
way, the new synthetic light curves preserve the noise properties of the original 
data. The fit procedure is repeated for each new light curve, eventually build- 
ing distributions in each of the fit parameters. We then estimate the la error 
bars from the analytic expression for the standard deviation of each distribution 
under the assumption that they are normally distributed. Each distribution is 
checked to ensure this assumption is valid. 

Figure 3 shows the bootstrap distributions for 10,000 realisations of the 
fit to our data. The phase is the most uncertain parameter. Amplitude and 
frequency distributions are normal in shape. Overplotted are Gaussians with the 
standard deviations derived from the distributions. The Gaussians are centred 
on the original fit parameters. It should be stressed that bootstrapping only 
estimates the uncertainties in parameters and does not refine them or assign 
significances to them. 

3. Discussion 

We determine that HD 134214 is monoperiodic and pulsating with a frequency 
of 2948.97 ± 0.55 ^iHz. This is consistent with the earlier photometric obser- 
vations of this star by Kreidl et al (1994). Those authors detected frequency 
changes which they interpreted as cyclic with a timescale of about 1 year. The 
frequency resolution of our short MOST run is insufficient to rigorously test 
this assertion. However, we have considered the implications of pulsational 
frequency variability in HD 134214. 

Heller & Kawaler (1988) have suggested that it may be possible to exploit 
the frequency variability of roAp stars to determine their evolutionary status. 
However, the one-year timescale of variability in the case of HD 134214 is too 
short to be associated with evolution. We have also modelled the amplitude 
of the frequency change over the main sequence life time of an A-type star. 
We estimate that a 20 year baseline of observations would show a frequency 
change on the order of 10~ 4 — 10~ 5 ^iHz. This is much smaller than the 
change of <~ 0.2 ^Hz reported by Kreidl et al (1994). Cameron et al. (in 
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preparation) have shown that taking magnetic perturbations to the oscillation 
frequencies into account does not improve the discrepancy between the observed 
and calculated frequency changes. 

Recent results by Kurtz et al. (2006) suggest that HD 134214 is oscillating 
spectroscopically in up to 6 modes. Our observations can rule out additional 
photometric oscillations at a 2a detection level of 0.36 mmag. Our estimated 
DFT noise level of ~ 0.12 mmag for 10 hours of observations with MOST is 
comparable to that obtained by Kreidl et al. (1994) based on about 56 hours 
of photometry from four observatories during about 4 months. The excitation 
and selection of pulsation modes in roAp stars is an open question. In the 
case of HD 134214, there are several avenues to be explored with respect to 
the additional modes seen spectroscopically: (1) Is it a case of radial velocities 
of certain elements and ionisation stages being more sensitive to degrees of 
higher 11 (2) If the sensitivity of the photometry can be improved sufficiently, 
will these modes also appear? (3) Is the broadband photometry just averaging 
over too wide an extent of the pulsating atmosphere? or (4) Might there be 
new physics in the upper atmospheres of (some) roAp stars to account for the 
differences observed? 

We estimate that even observations over only a few days with MOST would 
(in addition to improving our frequency resolution enough to investigate cyclic 
variability) reduce our noise levels by approximately half. Even this modest 
improvement in noise would be valuable as a test of the Kurtz et al. (2006) 
theory that photometric observations are not sensitive to the new type of upper 
atmosphere pulsational variability they report. 
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Figure 2: a) (top) The DFT of the HD 134214 data. Shown are the window function, 
the fit to the data (filled circle) and the bootstrap points (small dots). Also shown is 
the residual DFT after the fit; inverted for clarity, b) (bottom) The DFT of the data 
zoomed in on the region of the largest peak. Line styles are the same as in the top 
plot. 
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Figure 3: The bootstrap distributions for the fit parameters. Shown from top to 
bottom are the distributions for the phase <f> in radians, the amplitude in mmag, and 
the frequency u in /itHz. The labels above each panel give the mean n and the standard 
deviation a calculated for each distribution. Gaussians with width a and centered on 
the original fit parameters are plotted over each distribution. 



